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We consider a trapped atomic system in the presence of spatially varying laser fields. The laser- 
atom interaction generates a pseudospin degree of freedom (referred to simply as spin) and leads to 
an effective spin-orbit coupling for the fermions in the trap. Reflections of the fermions from the trap 
boundaries provide a physical mechanism for effective momentum relaxation and non-trivial spin 
dynamics due to the emergent spin-orbit coupling. We explicitly consider evolution of an initially 
spin-polarized Fermi gas in a two-dimensional harmonic trap and derive non-equilibrium behavior 
of the spin polarization. It shows periodic echoes with a frequency equal to the harmonic trapping 
frequency. Perturbations, such as an asymmetry of the trap, lead to the suppression of the spin 
echo amplitudes. We discuss a possible experimental setup to observe spin dynamics and provide 
numerical estimates of relevant parameters. 
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Ultracold atomic gases have proven to be an ideal 
test-ground for the experimental study of a variety of 
condensed matter phenomena A particularly in- 

teresting possibility is to realize spin-orbit (SO) in- 
teraction in cold atomic systems. In atomic gases, 
(pseudo)-spin represents a combination of different hy- 
pcrfine levels of atoms. It was proposed recently that mo- 
tion of atoms in position-dependent laser configurations 
may give rise to an effective non-Abelian gauge poten- 
tial 0, 0, 0, S HI ■ As argued in this Letter, a similar 
setup may lead to an effective SO interaction. In partic- 
ular, with proper laser configurations, one can engineer 
Rashba or linear Dresselhaus SO coupling terms as well 
as other types of couplings not accessible in solid state 
systems. The resulting spin dynamics in the adiabatic 
regime are remarkably simple, providing a robust experi- 
mental signature of the underlying effective SO coupling. 
Throughout this article, we refer to the pseudo-spin de- 
gree of freedom realized in cold- atom systems as simply 
spin. 

One of the main aspects of spin dynamics in bulk con- 
densed matter systems is the Dyakonov-Perel spin relax- 
ation [9j. This mechanism involves random elastic scat- 
tering of electrons off of impurities. These scatterings 
lead to spin relaxation, which is a result of spin preces- 
sion around a randomly oriented momentum-dependent 
axis. Another example is an electron in a quantum dot in 
the presence of SO coupling. In this case, the existence 
of discrete energy levels results in oscillatory time evo- 
lution of the spin polarization. The polarization decay 
and relaxation to an equilibrium state occur only in the 
presence of inelastic processes, such as electron-phonon 
interactions [l(J EH or nuclear hyperfine effects [H, [l]| . 
In the absence of such processes, the relaxation due to 
scatterings off of the quantum dot boundary is qualita- 
tively different from the bulk spin relaxation, as shown 
by the semiclassical analysis, e.g. of Chang et al. [3]. 



Atomic systems in the presence of spatially varying 
laser fields offer the new remarkable possibility of ob- 
serving strongly non-equilibrium quantum spin dynamics 
in a many-particle system without complications due to 
disorder that are always present in condensed matter sys- 
tems. In addition, by changing the geometry of the trap- 
ping potential one can access various dynamic regimes 
characterized by either regular or chaotic behavior of the 
spin polarization [lil flt| . In this Letter we concentrate 
on quantum spin dynamics of a Fermi gas confined in a 
two-dimensional harmonic trap and in the presence of a 
simple SO interaction term. We find that the spin po- 
larization may show two qualitatively different behaviors 
dependent on the strength of the SO coupling and on 
the number of particles. In the strong coupling regime, 
the initially fully polarized Fermi gas becomes completely 
unpolarized within a short time and remains unpolarized 
most of the time. In the weak coupling regime, the spin 
polarization never vanishes but oscillates with the pe- 
riod T — 2-k/lh, determined by the harmonic trapping 
frequency, u), with the oscillation amplitudes dependent 
on the SO coupling strength and the number of atoms. 
In both cases, spin polarization exhibits periodic echoes 
with the echo frequency equal to the harmonic trapping 
frequency. The echo amplitude and the period of spin po- 
larization oscillations are strongly modified in an asym- 
metric elliptic harmonic trap. 

Consider an ultracold Fermi gas confined in a quasi- 
two dimensional (xy-plane) harmonic trap. Along the z 
direction, the atomic dynamics is "frozen" by a hig h fre- 
quency optical trap [T3| or a deep optical lattice 18 1 . The 
SO coupling may be generated using the tripod scheme 
shown in Figflja), which is similar to the setup described 
in Ref. [Bj]. The states |1), |2), and |3) represent three de- 
generate hyperfine ground states (e.g., different Zeeman 
components of the hyperfine states of 6 Li atoms). These 
states are coupled to an excited state |0) by spatially 
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FIG. 1: Schematic representation of the experimental setup 
to create effective SO interaction, (a) Coupling between the 
lasers and the hyperfine states, (b) Energies of the degenerate 
dark states and of the bright states (\vi)). (c) Laser field 

configuration. 



varying laser fields L,o, with the corresponding Rabi fre- 
quencies ili, ^2, and SI3 . The single-particle Hamiltonian 



H = p 2 /2m + V ttap + Ha-t 



(1) 

where p is the momentum, Vtrap = J2jVj( r ) \j) (j\ 
represents the position-dependent trapping po- 
tential and Ti-a-i is the laser-atom interac- 
tion Hamiltonian, given by H a -i — A|0)(0| — 
[O 1 |0)(l|+O 2 |0)(2|+O 3 |0)(3|+h.c], where A is 
the detuning to the excited state |0). The Rabi fre- 
quencies can be parameterized as fii = fl sin 9 cos <f>e 1 , 
il 2 = fi sin 9 sin §e lS2 , and O3 = cos 9e lSi , with 

The diagonalization of 7i a _; 



yields four eigenstates: 



|1) 



two degenerate dark states \m) = sin</> e~ lSl3 

b e~ lS23 \2), \u 2 ) = cos6>cos0 e~' lSl3 |1) + 
)smd> e~ lS23 |2) — sin</>|3), and two non-degenerate 
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bright states \v ) « |0) and \v 3 ) w sin 9 cos <j) e lSl3 |1) 
sin sin e~ iS23 |2) - cos0|3), with = Si - Sj [Fig. 
[TJb)]. Note that the energy of the dark states is not 
modified by the laser field. We assume u> <C £1 <C A and 
oj <C ^ 2 /A, where w is the characteristic frequency of 
the trapping potential, and neglect contributions of or- 
der fi/A and smaller. This condition implies that the 
states from the subspace spanned by \u\) and \u2) are 
well separated in energy from the states \vq) and \va) , so 
that the coupling between dark and bright states is negli- 
gible (adiabatic approximation). In the present work we 
concentrate on a particular configuration of laser fields 
[see Fig[IJc)] characterized by Si = S2, S31 = S = mv s y 
, <fi — mv^x, and the constant angle 9 £ [0,7r]. The laser 
field Lio is generated by two laser beams propagating 
in the xy plane and intersecting at an angle, (. These 
two lasers form a standing wave along the n-direction 
and a plane wave along the y-direction. The laser field 
L20 has the same configuration except that the standing 
wave along the x-direction is shifted in phase by 7r/2. 
With such a setup, we have 4> = 2kix sin (C/2) [or mv^ = 
2ki sin (C/2)] and Si = S2 = fcij/cos (C/2), where ki = k 2 
are the wave-vectors for the laser beams. Finally, the 
laser field L30 is a plane wave propagating in the yz plane 



with an adjustable angle £ with respect to the y-axis, 
leading to S3 = k 3 y cos£ [mv s = fc 3 cos^ — ki cos (C/2))]. 

The effective low-energy Hamiltonian is obtained by 
projecting the Hamiltonian (TTJ) onto the subspace of the 
degenerate dark states \ui) ® |w 2 ) @| 



T~tu = 



P 

2m 



w(r) 



h + S a a z + Hso, 



(2) 



where I 2 represents the 2x2 unit matrix and &j with 
j £ {x, y, z} are the Pauli matrices. Here V\ — V 2 = w(r) 
and V3 = w(r) + S , where w(r) — muj 2 r 2 /2 (unless 
otherwise noted) is a symmetric harmonic potential and 
5 is a constant shift for the hyperfine state |3) that yields 
an effective Zeeman splitting 5q = sin 2 0[S— (v 2 +v 2 ^)/2]/2 
of the dark states and can be varied through additional 
state-dependent laser fields. [l9[ The last term in @ is 
the SO coupling 



7~Lso 



-v p x a y - vip y a z , 



(3) 



with vq — i>0cos# and vi — v s sin 2 9/2, where the 
pseudo-spin represents the internal degree of freedom as- 
sociated with the degenerate dark states \ui) and \u 2 )- 
We will call an atom in a \u\) (\u 2 )) state as a spin up 
(down) particle. The term Hso of the effective Hamilto- 
nian couples this (pseudo)-spin degree of freedom to the 
motion of the atom inside the trap. The "direction" of 
the spin does not have any real space significance and is 
solely determined by the choice of the basis in the dark 
state subspace. 

Atomic SO coupled systems open new possibilities to 
study strongly non-equilibrium spin dynamics, in con- 
trast to boundary linear response effects usually consid- 
ered in solid state systems with SO coupling 20]. Be- 
low, we study such a problem for the simplest Ising-like 
spin orbit coupling, which already leads to a non-trivial 
dynamics. This choice corresponds to a constant phase 
S = const; i.e., the plane- wave components of the three 
laser fields along y-direction share the same wave- 
vector. We define the spin polarization V{t) as the differ- 
ence between the occupation numbers of the dark states 
\u\) and \u 2 ). We assume that at t = 0, the system is 
fully polarized 'P(O) = 1, i.e. all the particles are in the 
dark state \ui). The corresponding single-particle quan- 
tum mechanics problem can be easily solved exactly. If 
</> a (r) is an eigenstate of the operator H = p 2 /2m + w(r) 
with an eigenvalue e a , then the eigenfunction of the SO 
coupled Hamiltonian H u — HqI 2 — vop x a y is 

iJ aX (r)=<j> a (r)e iXmV0X ^Q x ), (4) 

with A = ±1. For a harmonic trap, w(r) = mu> 2 r 2 /2, 
<^a(r) = l Pn x (x)ip rly (y) can be written as a product of 
the harmonic oscillator eigenfunctions tp n and the energy 
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FIG. 2: (Color online) Spin polarization as a function of time. 
Relaxation curves are plotted for TV = (N F + 1)(N F + 2)/2 
particles in a harmonic trap and in the presence of a (pseudo) 
SO interaction parameterized by a = (mv 2 /2u>) 1 ^ 2 . The 2tv/uj 
periodicity is due to the equal spacing between the harmonic 
oscillator levels. When aN^f 2 is large (black curve with tri- 
angles), the polarization is characterized by fast relaxation 
followed by periodic echoes. 



spectrum becomes e rix „ = to(n x + n y + 1). The eigen- 
functions ^ are degenerate with respect to A and have 
the eigenvalues E a — e a — mv 2 /2. 

Next, we define the spin polarization V(t) — 
jj($(t)\P z \&(t)) as an the average over the quantum 
many-body state <E>(i) obtained from the initial state 
$o by applying the time evolution operator U(t) = 
exp(— lYit). The operator P z can be expressed in terms 
of field operators ^(r) as P z — J d 2 r ^(r^^r) In 
the Heisenberg representation, the polarization becomes 
P z (t) — U~ x (t)P z U{t). We can obtain the explicit time 
dependence by expressing the field operators in terms of 
creation and annihilation operators associated with the 
single particle eigenstates (0]), <F(r) = J2 a ,X ^a\( r )^a\, 
and taking advantage of the simple time dependence of 
these operators, a a \(t) = a a \ e~ lEat . In order to eval- 
uate the action of the polarization operator on $0, it 
is convenient to introduce a set of operators b aa as- 
sociated with the fully polarized single particle states 

Vw( r ) = ^aWXtr, with x] = (1,0) and x\ = (0,1) 

and satisfying the relation b aa = J2/3\(' l Pacr\' , l'i3x)^a\- In 
terms of b-type operators, the initial state can be writ- 
ten as |$o) = 111 ^LtW wnere |0) represents the vacuum 
and the prime signifies that the product is constrained 
by the condition E a < Ep, where Ep = u(Np + 1) is the 
"Fermi" energy. After recasting P z (t) in terms of b-type 
operators, the polarization becomes 



(E^<E F ) 



V(t) 



1 

2N 5I„ 



(5) 



<r,a,f3 7 



i(e Q -e/3)t 



where ^ = mv^x and a = ±1. Eq. ^ expresses the 
polarization of a non-interacting many-body system with 



FIG. 3: (Color online) The time independent contribution to 
the polarization as a function of the effective coupling param- 
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SO interaction in terms of matrix elements of the single 
particle states 4> a describing the motion of an atom in 
the confining potential. The effects of the SO coupling 
are contained in the parameter £ and the finite number 
of particles enters through the constraint on the 7 sum. 
Equation ([5]) contains matrix elements (<p„|e~ 2lcr ^ \f m ) 
of the harmonic eigenfunctions tp n (%) with energies e„ = 
uj(n+ 1/2). The summations over the quantum numbers 
are performed numerically. The time dependence of the 
pseudo-spin polarization for N particles that arc initially 
in the dark state and occupy the first Np energy 
levels is shown in Fig. [51 We parameterize the strength 
of the SO coupling by a — (mi;g/2?iw) 1 / 2 . Notice that 
all the curves in Fig. [2] display a periodic structure. 

Spin relaxation depends strongly on both the strength 
of the SO interaction (i.e. on a) and the number of parti- 
cles (i.e. on Np). In the weak coupling limit, a 2 Np <C 1, 
the polarization deviates slightly from unity and we have 



T{t) w 1 - 4a 2 



+ 1 ) [1 - cos(wt)]. 



(6) 



At strong coupling, a 2 Np ^> 1, the polarization is char- 
acterized by fast spin relaxation to zero and periodic spin 
echoes that restore V to its initial value. Performing a 
power expansion of the polarization at small times we get 



(7) 



where c„ are numerical coefficients. From this we obtain 
the time scale of the fast relaxation, At = l/(ujayNp), 
which is also the characteristic width of the echo peaks. 

The relaxation to zero polarization in the strong cou- 
pling limit is a non-trivial dynamical effect. In Eq. ((SJ, 



we can separate explicitly the time independent contribu- 
tions coming from the states with e a = ep and write the 
polarization as V{t) = P CO nst + AP(t). For large but fi- 
nite values of the parameter ct\/Np the time-independent 
contribution P CO nst does not vanish. Instead, after the 
initial relaxation, the time-dependent term approaches 
a constant value, AP(i) — * — P CO nst, leading to a van- 
ishing polarization. The time independent contribution 
Pconst is shown in Fig. [3] for different sets of parameters 
(a,Np). Notice that in the large coupling limit P CO nst 
depends on a single scaling parameter, oi\J Np. In fact, 
as it is suggested by Eq. (|Tj) , the polarization itself scales 
with a^Np in the strong coupling limit (inset of Fig. [3]). 

The periodic time dependence of the polarization is a 
consequence of the equidistant energy spectrum. What 
would be the result of altering this structure? We con- 
sider an elliptical trap described by the potential w(r) = 
m\jj\(x cosx — J/sinx) 2 + sin x + y cosx) 2 ]/2, where 
the angle x describes the orientation of the symmetry 
axes of the trap potential relative to the (x, y) directions 
defined by the laser field. If \ — (x — 7r /2), P(i) is the 
same as for an isotropic system with a trap frequency lui 
(1^2)- If the orientation of the symmetry axes deviates 
from the (x,y) directions, V(t) is still periodic if ^2/^1 
is a rational number and non-periodic otherwise. In the 
former case with 0J2/0J1 = p/q, the period of oscillations 
is T = p(2-k/lu 2 ) = q(2n/uj 1 ) (see, Fig. g]). Notice that 
in the weak coupling regime (upper panel in Fig. 2]) for 
small deviations from the symmetry axes, spin relaxation 
is reminiscent of the periodic structure of the symmet- 
ric configuration. However, the residual echoes become 
weaker and eventually vanish in the strong coupling limit 
(lower panel in Fig. [3J). 

Finally, we briefly discuss how to observe spin relax- 
ation and echoes in ultracold Fermi gases through the 
time of flight measurements. Initially the atoms are pre- 
pared in the state To transfer them to the spin up 
polarized state one applies a Raman pulse between 
the states |1) and |2) with a spatially dependent Rabi 
frequency that matches the spatial variation of the \ui) 
state. Then, the laser fields Lio are turned on and the 
Fermi gas experiences the effective SO coupling. After 
a certain time i, one turns off the laser fields Liq and 
applies a reversal Raman pulse to transfer the atoms 
from \ux) back to |1). Finally, a time of flight measure- 
ment gives the number of atoms in the |1) state (i.e. in 
at time t) and thus determines the spin polarization 
P(i). A possible choice of parameters is 10 — 2tt x 10Hz, 
A = 2?r x 10 2 GHz , and Q = 2tt x 10 7 Pz, satisfying 
lu < CI < A and to < f2 2 /A. The strength of the SO 

coupling vq = v^cosO — ^1 sin cos# can be var- 
ied from to ^Ml cos 9 by adjusting the angle C between 
the laser beams. Consequently, a — (rav^^UuS) 1 ! 2 — 
2 (ujr/uj) 1 / 2 sin cos# (where ui r is the photon recoil 




FIG. 4: (Color online) Spin polarization as a function of time 
in an elliptic harmonic trap (with lj2/<^i = 23/33). Relax- 
ation curves are plotted for different laser orientations pa- 
rameterized by the angle \ ( see text). The upper and lower 
panels correspond to the SO couplings a = 0.2 and a = 0.6, 
respectively. For an arbitrary orientation of the lasers, V(t) 
has a period 23(27r/tiJ2) = 33(27r/o;i). The residual echoes 
become smaller and vanish in the strong coupling limit. 



frequency for Lio), can vary in a range from to 100 
for 6 Li atoms. Therefore, both the strong and weak SO 
coupling limits are accessible in experiment. 

C.Z. was supported by ARO-DTO, ARO-LPS, and 
LPS-NSA. 
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